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Abstract 

We calculate the Next-to-Leading Order (NLO) QCD corrections to Wbb production includ- 
ing full bottom-quark mass effects. We study the impact of NLO QCD corrections on the total 
cross section and invariant mass distribution of the bottom-quark jet pair at the Fermilab Teva- 
tron pp collider. We perform a detailed comparison with a calculation that considers massless 
bottom quarks. We find that neglecting bottom-quark mass effects overestimates the NLO to- 
tal cross-section for Wbb production at the Tevatron by about 8% independent of the choice of 
renormalization and factorization scale. 
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I. INTRODUCTION 



The associated production of a boson with a bb pair plays a critical role at the Fermilab 
Tevatron pp collider, since it accounts for one of the most important background processes 
to both the associated production of a Higgs boson with a W boson, pp — > HW (with 
H ^bb) i Q, S, y, Q, and single-top production, pp tb, ib (with t{t) Wb{b)) ^,0,3. 
These two processes are of extreme relevance to the physics program of the Tevatron: they 
both test fundamental predictions of the Standard Model (SM), i.e. the existence of a Higgs 
boson and the structure of the Wtb vertex, and at the same time constitute a window to 
new physics. The cross section for pp HW has been calculated including up to Next-to- 
Next-to-Leading (NNLO) QCD corrections and 0{a) electroweak corrections 

while single-top production has been calculated at Next-to-Leading (NLO) in QCD 
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21[, and at one-loop of electroweak (SM and MSSM) corrections 



IS |2^. 



The production of a Higgs boson in association with an electroweak gauge boson, pp — > 
HV {V = Z, W) with H bb, is the most sensitive production channel of a SM Higgs 
boson at the Tevatron for a Higgs boson lighter than about 140 GeV. A relatively light SM 
Higgs boson is preferred by electroweak precision data, Mh = 89ttl GeV at 68% confidence 
level 231 ^. The Tevatron with an integrated luminosity of 6 fb~^ will be able to exclude 

n 

a Higgs boson with 115 GeV < Mh < 180 GeV at 95% confidence level j2i], which will 
provide important guidance for the search strategy at the CERN Large Hadron Collider. 
Thus, to fully exploit the Tevatron's potential to detect the SM Higgs boson or to impose 
limits on its mass, it is crucial that the dominant background processes are under good 
theoretical control. 

In the present experimental analyses^ the effects of NLO QCD corrections on the total 
cross-section and the dijet invariant mass distribution of the Wbb background process have 
been taken into account by using the MCFM package In MCFM, the NLO QCD 

predictions of both total and differential cross-sections for the qq' —>■ Wbb production process 
have been calculated in the zero bottom-quark mass (m^ = 0) approximation [2^ 23|. From a 
study of the Leading Order (LO) cross-section, finite bottom-quark mass effects are expected 



^ For an update see the LEPEWWG website at |http: / /lepewwg.web.cern.ch/LEPEWWG 

^ For updated results, see the CDF and D0 websites at www-cdf.fnal.gov/physics/exotic/exotic.html and 

www-dO.fnal.gov/Run2Physics/WWW/results/higgs.htm. 
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to affect both total and differential Wbb cross-sections mostly in the region of small 66-pair 



invariant masses 27[ . Given the variety of experimental analyses involved both in the search 
for HW associated production and single-top production, it is important to precisely assess 
the impact of a finite bottom-quark mass over the entire kinematical reach of the process, 
including complete NLO QCD corrections. 

In this paper we compute the NLO QCD corrections to qq' — Wbb, including full bottom- 
quark mass effects. Using the MCFM package 2^, we compare our results with the cor- 
responding results obtained in the rrih = limit. Numerical results are presented for the 
total cross-section and the invariant mass distribution of the bb jet pair, for the Tevatron pp 
collider, including kinematic cuts and a jet-finding algorithm. In particular, we apply the 

jet algorithm and require two tagged 6-jets in the final state. 

The outline of the paper is as follows. In Section |n] we briefly discuss the technical details 
of our calculation, while we present numerical results and a discussion of the bottom-quark 
mass effects in Section IIIII Section HXl contains our conclusions. 



II. CALCULATION 

The NLO QCD corrections to qq' —>■ bbW consist of both one-loop virtual corrections to 
the tree level processes depicted in Fig. ^ and one-parton real radiation from both the initial 
and final state quarks, i.e. qq' — » bbW + g. At the same order, the qg{qg) — ^ bbW + q{q) 
process also needs to be included. 
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FIG. 1: Tree level Feynman diagrams for 



bbW. 



The 0{as) virtual corrections consist of self-energy, vertex, box, and pentagon diagrams 
with several massive propagators, since we take nib ^ 0. They contain both UV and IR 
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FIG. 2: Dependence of the total NLO QCD cross-section on the 6s PSS parameter, when 5c is 
fixed at 5c = 10~^. In the upper window we illustrate separately the cutoff dependence of the 
soft and hard-collinear part (2 3, red dashed curve) and of the hard non-collinear part (2 — s- 4, 
blue dotted curve) of the real corrections to the total cross-section. The 2 — > 3 curve also includes 
those parts of the 2^3 NLO cross-section that do not depend on 5c and 5s, i.e. the tree level 
and one-loop virtual contributions. The sum of all the contributions corresponds to the black solid 
line. The lower window shows a blow-up of the black solid line in the upper plot, to illustrate 
the stability of the result. The error bars indicate the statistical uncertainty of the Monte Carlo 
integration. 



singularities which need to be computed analytically. For box an d p entagon diagrams we 
use techniques similar to the ones explained in detail in Refs. js^^Sl- dimensional 
regularization with d = 4— 2e, and extract both UV and IR divergences as poles in e. The UV 
singularities are cancelled by introducing a suitable series of counterterms. We renormalize 
the wave functions of the external quark fields in the on-shell scheme, and the strong coupling 
constant in the MS-scheme, decoupling the top quark. At this order in QCD the weak 
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FIG. 3: Dependence of the total NLO QCD cross-section on the 5c PSS parameter, when 5s is 
fixed at 5s = 10^^. In the upper window we illustrate separately the cutoff dependence of the 
soft and hard-collinear part (2 3, red dashed curve) and of the hard non-collinear part (2 — s- 4, 
blue dotted curve) of the real corrections to the total cross-section. The 2 — > 3 curve also includes 
those parts of the 2^3 NLO cross-section that do not depend on 5c and 5s-, i.e. the tree level 
and one-loop virtual contributions. The sum of all the contributions corresponds to the black solid 
line. The lower window shows a blow-up of the black solid line in the upper plot, to illustrate 
the stability of the result. The error bars indicate the statistical uncertainty of the Monte Carlo 
integration. 



vertex renormalization consists only of the external quark wave-function renormalization. 
Self-energy, vertex, box, and pentagon diagrams contain IR divergences that combine and 
cancel against the analogous divergences in the real emission 0{as) corrections, and in the 
renormalized Parton Distribution Functions (PDF). 

We compute the real emission 0{as) corrections using the Phase Space Slicing (PSS) 
method with two cutoffs: 5s for the soft singularities, and 5c for the hard-collinear singu- 
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30|. The independence of the final result on the arbitrary values of Sg and 



6c has been checked over a large range of values for both parameters and is illustrated in 
Figs. 121 and El The numerical results in Section UTTl have been obtained using 6s = 10^^ and 

4 = 10-^ 

In our calculation we treat 75 according to the naive dimensional regularization approach, 
i.e. we enforce that 75 anticommutes with all other 7 matrices in c? = 4 — 2e dimensions. This 
is known to give origin to inconsistencies when at the same time the ci-dimensional trace of 
four 7 matrices and one 75 is forced to be non-zero (as in d = 4, where Tr(7'^7'^7''7'^75) = 
Aie^'^''") [2^. In our calculation both UV and IR divergences are handled in such a way 
that we never have to enforce simultaneously these two properties of the Dirac algebra in d 
dimensions. For instance, the UV divergences are extracted and cancelled at the amplitude 
level, after which the d ^ A limit is taken and the renormalized amplitude is squared using 
d = A. Thus, all fermion traces appearing at this point are computed in four dimensions 
and therefore have no ambiguities. 

Both virtual and real corrections have been checked by independent calculations that have 
used FORM Q, TRACER H, 

the FF package Q, and MAPLE_^he,2 ^ 4 amplitudes 
for the real corrections have been double checked using Madgraph 



III. NUMERICAL RESULTS 

In this paper we present results for Whh production at the Tevatron, including NLO 
QCD corrections, and using a non-zero bottom-quark mass, fixed at m6=4.62 GeV. The W 
boson is considered on-shell and its mass is taken to be My/ = 80.41 GeV. The mass of 
the top quark, entering in virtual corrections, is set to = 174 GeV. The LO results use 
the 1-loop evolution of as and the CTEQ6L set of PDF j^^, while the NLO results use 
the 2-loop evolution of a, and the CTEQ6M set of PDF, with a^^^{Mz) = 0.118. The 
W boson coupling to quarks is proportional to the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. We take Vud = = 0.975 and Vus = Vcd = 0.222, while we neglect the 
contribution of the third generation, since it is suppressed either by the initial state quark 



densities or by the corresponding CKM matrix e 



We implement the kx jet algorithm 



41 



ements. 



42] with a pseudo-cone size i? = 0.7 and we 



recombine the parton momenta within a jet using the so called covariant ii^-scheme j4Q|. We 
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checked that our implementation of the /ct jet algorithm coincides with the one in MCFM. 
We require all events to have a bb jet pair in the final state, with a transverse momentum 
larger than 15 GeV (j9^^ > 15 GeV) and a pseudorapidity that satisfies \ri^'^\ < 2. We impose 
the same px and |?7| cuts also on the extra jet that may arise due to hard non-coUinear real 
emission of a parton, i.e. in the processes Wbb + g or Wbb + q{q). This hard non-collinear 
extra parton is treated either inclusively or exclusively, following the definition of inclusive 
and exclusive as implemented in the MCFM code j25|. In the inclusive case we include 
both two- and three-jet events, while in the exclusive case we require exactly two jets in the 
event. Two-jet events consist of a bottom-quark jet pair that may also include a final-state 
light parton (gluon or quark) due to the applied recombination procedure. Results in the 
massless bottom-quark approximation have been obtained using the MCFM code 

5 I 1 < 1 < < ^ 
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FIG. 4: Dependence of the LO (black solid band), NLO inclusive (blue dashed band), and NLO 
exclusive (red dotted band) total cross-sections on the renormalization/factorization scales, includ- 
ing full bottom-quark mass effects. The bands are obtained by varying both ^u/j and fip between 
/io/2 and 4^o (with = n^-b + -^vk/2). 



In Figs. Elini we illustrate the renormalization and factorization scale dependence of the 
LO and NLO total cross-sections, both in the inclusive and exclusive case. Fig. |3] shows 
the overall scale dependence of both LO, NLO inclusive and NLO exclusive total cross- 
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FIG. 5: Dependence of the LO and NLO inclusive total cross-section on the renormahza- 
tion/factorization scale, when /i/j = hf- The left hand side plot compares both LO and NLO 
total cross-sections for the case in which the bottom quark is treated as massless (MCFM) or 
massive (our calculation). The right hand side plot shows separately, for the massive case only, the 
scale dependence of the qq' and qg + qg contributions, as well as their sum. 

sections, when both and hf are varied independently between /io/2 and 4/io (with /iq = 
mb + Miy/2), including full bottom-quark mass effects. We notice that the NLO cross- 
sections have a reduced scale dependence over most of the range of scales shown, and the 
exclusive NLO cross-section is more stable than the inclusive one especially at low scales. 
This is consistent with the fact that the inclusive NLO cross-section integrates over the 
entire phase space of the qg{qg) bbW + q{q) channels that are evaluated with NLO as 
and NLO PDF, but are actually tree-level processes and retain therefore a strong scale 
dependence. In the exclusive case only the 2-^3 collinear kinematic of these processes is 
retained, since 3-jets events are discarded, and this makes the overall renormalization and 
factorization scale dependence milder. To better illustrate this point, we show in the right 
hand side plots of Figs. El and IHl the mu-dependence of the total cross section and of the 
partial cross-sections corresponding to the qq' and the qg + qg initiated channels separately, 
for /i/j = fiF, both for the inclusive and for the exclusive case. It is clear that the low scale 
instability of the inclusive cross-section is entirely driven by the qg + qg contribution. In 
the left hand side plots of Figs. El and El we also compare the scale dependence of our results 
to the scale dependence of the corresponding results obtained with ruh = (using MCFM), 
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FIG. 6: Dependence of the LO and NLO exclusive total cross-section on the renormahza- 
tion/factorization scale, when /i/j = fJ-F- The left hand side plot compares both LO and NLO 
total cross-sections for the case in which the bottom quark is treated as massless (MCFM) or 
massive (our calculation). The right hand side plot shows separately, for the massive case only, the 
scale dependence of the qq' and qg + qg contributions, as well as their sum. 



both at LO and at NLO. Using a non-zero value of is not expected to have any impact 
on the scale dependence of the result^ and, indeed, the scale dependence of the LO and NLO 
pair of curves is very similar, with a shift due to the bottom-quark mass effects. 

While the LO cross-section still has a 40% uncertainty due to scale dependence, this un- 
certainty is reduced at NLO to about 20% for the inclusive and to about 10% for the exclusive 
cross-section respectively. The uncertainties have been estimated as the positive/negative 
deviation with respect to the mid-point of the bands plotted in Fig. EJ where each band range 
is defined by the minimum and maximum value in the band. We notice incidentally that 
the difference due to finite bottom-quark mass effects is less significant than the theoretical 
uncertainty due to the residual scale dependence in the inclusive case, but is comparable in 
size in the exclusive case. Indeed, the finite bottom-quark mass effects amount to about 8% 
in both inclusive and exclusive cases. 

In Fig. 13 we show the rescaled difference between the total cross-sections obtained from 

^ Note that we always use rrif, = 4.62 GeV in the determination of the scales in terms of /ip ~ rni, + Mi^'/2 
even in the results obtained with mj, — 0. 
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FIG. 7: Dependence on the renormalization/factorization scale of the rescaled difference between 
our NLO calculation (with ^ 0) of the total cross-section and MCFM (with mf, = 0) for the 
inclusive and exclusive cases (with hr = ij,f)- The error bars indicate the statistical uncertainty of 
the Monte Carlo integration. 



our calculation (with ^ 0) and with MCFM (with m;, = 0) defined as follows: 

A. . .-o„„, ^ 0) - ^ 0) . 

As can be seen, within the statistical errors of the Monte Carlo integration, the finite bottom- 
quark mass effects on the total cross-sections at NLO are well described by the corresponding 
effects at LO. 

Finally, in Figs. ISHT^we study the distribution da/dmi,i, where m^^ is the invariant mass 
of the bb jet pair. The impact of NLO QCD corrections on this distribution is illustrated in 
Figs. IHl and ini for the inclusive and exclusive case respectively. We see that the NLO QCD 
corrections affects the cross section quite substantially in particular for low values of rrii^i. In 
each figure the right hand side plot gives the ratio of the NLO and LO distributions, providing 
a sort of i^-factor bin by bin. Figs. ITUl and ITD compare the NLO rfcr/rfm^^ distributions 
obtained from the massive and massless bottom-quark calculations. The results with m^, = 
have been obtained using MCFM. As expected, most of the difference between the massless 
and massive bottom-quark cross-sections is coming from the region of low invariant mass 
m^ft, both at LO and at NLO, where the cross-sections for nib are consistently below the 
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FIG. 8: The inclusive distribution da/dm^j^ in LO and NLO QCD. The right hand side plot shows 
the ratio of the LO and NLO distributions. 

ones with mf, = 0. For completeness, we also show in Fig.^lthe comparison between massive 
{nib 7^ 0) and massless {mb = 0) calculations at LO in QCD. The LO m^g distribution for 
massive bottom-quarks has been obtained both from our calculation and from MCFM, which 
implements the rrib option at tree level, and both results have been found in perfect 
agreement. As can be seen by comparing Figs. ITUlfTTl and Fig. the impact of a non-zero 
bottom-quark mass is almost not affected by including NLO QCD corrections. To illustrate 
this in more detail we show in Fig. El the rescaled difference between the rribb distributions 
obtained with our NLO calculation (with rrib 7^ 0) and with MCFM (with nib = 0) defined 
as follows: 

, da da'^'^o da^'^o , , rf(T^«(mb ^ 0) 

^d^b ^ ^^"^ ^ " ^^"^ ^ da4mb - 0) • 
Apart from small deviations in the rubb region below about 100 GeV, the finite bottom-quark 

mass effects at NLO are well described by the LO calculations. 

IV. CONCLUSIONS 

We have calculated the NLO QCD corrections to qq' Wbb production including full 
bottom-quark mass effects. We have presented numerical results for the total cross-section 
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and the invariant mass distribution of the bottom-quark jet pair (m^j) at the Tevatron for 
both massless and massive bottom quarks. We apply the kx jet algorithm, require two h- 
tagged jets, and impose kinematical cuts that are inspired by the D0 and CDF searches for 
the SM Higgs boson in WH production. The bottom-quark mass effects amount to about 
8% of the total NLO QCD cross-section and can impact the shape of the m^^ distributions, 
in particular in regions of low mf,^. This is relevant to SM Higgs searches in WH associated 
production and to searches for single-top production. 
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FIG. 13: The distribution of the rescaled difference between our NLO calculation (with rrib 
and MCFM (with mf, = 0) for the inclusive (upper plot) and exclusive case (lower plot). 
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